Nature of insulating state in NaV205 above charge-ordering transition: a cluster 

DMFT study 
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The nature of insulating state driven by electronic correlations in the quarter-filled ladder com- 
pound a'NaV205 is investigated within a cluster dynamical mean-field approach. An extended 
Hubbard model with first-principle tight-binding parameters have been used. It is shown that the 
insulating state in the charge-disordered phase of this compound is formed due to the transfer of 
spectral density and dynamical charge fiuctuations where for the latter, the role of inter-site Coulomb 
interaction is found to be of crucial importance. 



The ladder compound NaV205 has been a subject of 
great interest since last five years ]l|-|l^. It exhibits a 
remarkable phase transition at Tc =34 K, now identified 
as charge-ordering of zigzag type |g|,||-||]. Both charge- 
ordered phase and charge-disordered one are insulating 
with an energy gap of the order of 0.8-1 eV H. The 
presence of the gap in the ordered phase is not surpris- 
ing and it was reproduced successfully, for example, in 
the recent local density approximation(LDA)-fU calcu- 
lations H. The properties of the disordered phase are 
much more difficult to understand. In contrast to the 
iso-structural ladder compound, CaV205, NaV205 has 
quarter-filled band |0] rather than half-filled one and 
cannot be considered as a standard Mott insulator pTl. 
It has been proposed in Ref. H that a large value of the 
transverse hopping parameter in the ladder which splits 
the band into sub-bands of the bonding and anti-bonding 
states could be responsible for the Mott insulator be- 
haviour in NaV2 05 in the presence of strong Coulomb 
interactions. However it is not known whether this mech- 
anism is adequate for the realistic values of the parame- 
ters characterising the single-particle electronic structure 
and the electron-electron correlations. 

In this letter we investigate the correlation effects in 
NaV205 taking into account non-local dynamical charge 
fluctuations on the rung for this two-leg ladder com- 
pound. This allows us to understand the origin of the 
insulating states above Tc and to estimate the relative 
importance of various physical mechanisms responsible 
for the gap formation. 

The crystal structure of NaV205 projected in the xy 
plane is schematically shown in Fig.l. The results of the 
X-ray |q| , NMR m , and optical 0] experiments as well 
as the Hartree-Fock calculations ^ support the zigzag 
charge-ordering state for low temperatures (see Fig.l). 
In this state one has approximately one d electron per 
rung of the vanadium ladder. We start with LDA-I-U ||l^] 



calculations of the ordered states but in contrast to the 
previous work H we considered several different types 
of the charge ordering. This gives us an opportunity 
to estimate the on-site and inter-site Coulomb interac- 
tions U and V respectively which in turn were used to 
parametrise the model Hamiltonians to be used for the 
calculations taking into account the dynamical correla- 
tion effects. 




FIG. 1. Schematic representation of the crystal struc- 
ture of the vanadium layers in NaV2 05 and the hopping 
matrix elements. The vanadium ions are denoted by filled 
circles. The ellipse shows the cluster which plays the 
role of an effective impurity in the DMFT calculations. 
A zig-zag charge ordering of the V*^ and V^""" ions, ob- 
tained from our LDA-(-U calculations as a ground state is 
shown on the left ladder. Bold arrows are the translation 
vectors. 



By mapping the results of the LDA+U calculations for 
different types of charge ordering on the results of model 
calculations with on-site (U) and inter-site (V) Coulomb 
interaction parameters, we obtained the following values: 
U=2.8 eV and V=0.17 eV. 

It is natural to assume that the tendency to keep the 
number of d electrons per rung close to unity also takes 
place above the transition temperature leading to strong 
short-range order and well-developed dynamical charge 
fluctuations. This is confirmed by the temperature de- 
pendences of the spin gap and the entropy measurements 
pO| . Usual LDA as well as the mean-field theories like 
Hartree-Fock or LDA-I-U methods are insufficient to take 
into account these essential many-body processes. The 
simplest reliable way to consider such short-ranged corre- 
lation effects is the use of the dynamical mean-field the- 
ory (DMFT) |l^ which can be combined with realistic 
LDA band structure calculations (LDA+DMFT) @-0. 
The DMFT maps the initial many-body problem for a 
crystal onto a self-consistent quantum-impurity problem. 
To consider the phenomena such as charge ordering or 
fluctuations, the inter-site correlations are of crucial im- 
portance, therefore we need to use a cluster generalisa- 
tion of the DMFT method |l7|-|l|]. The most reasonable 
choice of the cluster in our case is a pair of vanadium 
atom at the rung (see Fig.l). The unique geometry of 
the ladder compounds makes the choice of proper clus- 
ter simpler and the free cluster consideration is the most 
natural choice p^Jl9||. 
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FIG. 2. Bare densities of states (DOS). Solid lines cor- 
respond to LDA DOS. The dashed and long dashed lines 
represent the first-principle tight-binding parameterisa- 
tions without and with td hopping respectively. 



We use the extended Hubbard model for two-leg lad- 
der: 

ij i *<i,cr 



where Uj-s are the effective hoppings, Ui and Vij are local 
and inter-site Coulomb interactions respectively, riia- = 
cf^Cia- The proper choice of hopping parameters is not 
simple, and the most widely used set ( t± = 0.38 eV, 
t|[ = 0.18 eV, ii = 0.012 eV, ia = 0.03 eV ) was obtained 
by the fitting to the LDA bands 12]. Recently a rigor- 
ous procedure of massive downfolding of LDA bands to a 
few-band description and the subsequent Fourier trans- 
formation of the resulting Hamiltonian from the recip- 
rocal to direct space to extract the single-electron pa- 
rameters has been developed within the framework of 
linear-muffin-tin-orbital (LMTO) description pOl. This 
method applied to NaV205 gave the following set of hop- 
ping parameters: t_L=0.398 eV, t|| =0.084 eV, ti=0.025 
eV, t2=0.022 eV which is rather close to the standard 
one presented above, but in addition the diagonal (Fig. 
1) hopping parameter td=0.083 eV is appreciable and 
was not considered before. We have found that these di- 
agonal hopping processes are very important. As can be 
see in Fig. 2, the inclusion of diagonal hopping td in the 
single-electron part of the model Hamiltonian provides a 
much better agreement of the bare DOS with the LDA 
DOS. It also results in the insulating state of the charge- 
disordered systems for the realistic values of U and V 
obtained from LDA-f U calculations, while the other TB 
parametrisation, excluding the diagonal hopping, results 
into the metallic state for the same values of U and V 
(see Fig. 3). 

The crystal Green function matrix in LDA+DMFT ap- 
proach can be written as 

G (k,iw) = [{iuj + /i) * 1 - h (k) - S iiu)r^ 

where hap (k) is the effective hopping matrix, S^^ (icu) 
is the self-energy matrix of the two-site super-cell dimen- 
sion which is assumed to be local, i.e. k-independent, 
and /i is the chemical potential. 

In the cluster version of the DMFT scheme |l^,|l3| , one 
can write the matrix equation for a bath Green function 
matrix Q which describe an effective interactions with the 
rest of the crystal: 



g-^ {iu) = G^ (iuj) + S (iu) 



(1) 



where the local cluster Green function matrix is equal to 
Gaf3 (iw) = Y^ Gaf3 (k,iu;) , and the summation runs over 

k 

the Brillouin zone of the lattice. 

We used the exact diagonalisation (ED) scheme to 
solve the cluster DMFT problem. In this case the lat- 
tice Hamiltonian is mapped onto finite cluster impurity 
model: 

Htmp= ^ T,jcl^Cj„ + ^ U.,ni-^nii + ^ VijU^^n.^a (2) 
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k,i,j,(T 



where T^ is the hopping parameters inside the cluster 
(for our two site rung-cluster this is only t±), Eija{k), 
and Tij„{k) are effective energies and hybridisation ma- 
trix for finite-chain bath orbitals k = 1, ...,nh. 

For the iterative solution of the effective impurity 
model Eq. (||) we use the Lanczos version of ED method 
p3[ . The orbital energy matrix for the conduction band 
Eija{k), and the corresponding hybridisation elements 
rjjo-(fc), are the effective parameters which reproduce the 
bath Green function: 



(3) 

In the paramagnetic case we can transform Q, T, r^, 
and Efe matrices of the dimension 2x2 for our clus- 
ter to the diagonal bonding- antibonding basis A = {b,a}. 
In this case Tx = {—t±,t±} and Q\ = {gb,Ga}, where 
Gb,a = {Gil ± 5i2)/2. We used 10 bath orbitals for our 
two-site cluster. 

The parameters {EkXaT^kXa} are now fitted to repro- 
duce the bath Green function Gxaii'jJv) for bonding and 
anti-bonding states independently. The next step is the 
solution of the cluster-impurity problem (Eq. (^) to 
get the cluster self-energies Tnj„{iLOi,) which are required 
for the next DMFT- iterations. After solving the effec- 
tive cluster problem using Lanczos algorithm, the local 
Green's function matrix Gija{iuJu) was determined and 
the cluster self-energy T,ija{iiOv) was obtained within the 
cluster- impurity scheme (Eq. d3)). 

The phase digram obtained in our DMFT calculations 
is shown in Fig. 3. One can see that for large enough U 
and realistic values of V the disordered state turns out 
to be insulating. This arises from the two physically dif- 
ferent mechanisms. The first one is the spectral density 
transfer which can effectively change the quarter-filled 
system into the half- filled one. For large enough U the 
energy band splits into Hubbard sub-bands with the av- 
erage spectral weight 1/2 for each of them and therefore 
the system appears to be insulator for half-filling instead 
of complete filling |£1|] . In this large- U limit the Hubbard 
model reduces to so-called t-J-V model and can explain 
insulating properties and optical spectra of NaV205 p2[ . 
Our calculations do show a formation of the lower and 
upper Hubbard bands which leads to the spectral density 
transfer and moves the Fermi energy into the pseudo-gap 
between bonding and anti-bonding states in the lower 
Hubbard band as it was proposed qualitatively in Rcf. 
H. However, the gap which can be obtained due to this 
mechanism is very small. To increase its value (consistent 
with the experiment a) the inclusion of the inter-site 
Coulomb repulsion (V) appears to be important (Fig. 
4). Our calculations show that the broad enough gap 
arises with the increase of V suddenly, as a result of a 
first-order phase transition. With V=0 the insulator is 



stable only for U « 4 eV and above, while already small 
value of V=0.1 eV decreases the critical U value below 3 
eV. The resuUs of LDA-I-U calculations (U=2.8 eV and 
V=0.17 eV) gave a point in the phase diagram (the cross 
on Fig. 3) which is well above the metal-insulator transi- 
tion line. It is important to stress that if one carries out 
the calculation in usual single-site DMFT calculation in- 
stead of the cluster DMFT an adequate description of 
the electronic structure is not obtained for the same val- 
ues of the parameters, and the insulating states appears 
only for U > 12 eV! This demonstrates the crucial im- 
portance of the charge fluctuations on the rung for the 
formation of insulating state in NaV205 with the realistic 
U- values. We also analysed the structure of the ground 
state in the effective cluster model (Eq.(||)) and found 
that for our LDA-I-U Coulomb parameters about 70% of 
the ground state eigenfunction corresponds to the config- 
urations with one electron per rung and 20% related to 
empty rung states. This means that the configurations 
with the empty rungs in NaV205 has also appreciable 
weight. 
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FIG. 3. Calculated phase diagram for the model with 
hopping parameter from Ref. H (dashed line) and the pa- 
rameters obtained from our TB-LMTO calculation (solid 
line). The lines demark the metallic and insulating 
phases in two calculations. The cross corresponds the 
values of U and V obtained from our LDA-fU calcula- 
tion. 



We conclude that the insulating state of NaV205 above 
T =34 X is characterised by strong dynamical charge 
fluctuations. In this respect the situation is close to 
the half-filled Hubbard chain or an pseudospin-liquid like 
phase: the tendency for the formation of the state with 
one electron per rung is similar to the singlet state for- 
mation in strongly frustrated spin systems and can be 
described in the limit of large U and V by anisotropic 
Heisenberg-like model p^,0,p|,p4l. It is important to 



stress that experimentally the spin gap exists also above 
Tc pX(|. However, this description is only qualitative. 
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FIG. 4. The renormalised densities of state (only the 
lower Hubbard bands) obtained in our DMFT calculation 
with parameter U=2.8 eV and V^=0.17 eV(solid line); 
17=2.8 eV and l/=0.5 eV(dashed line). 

Our calculations demonstrate that this compound is 
situated near the boundary between metal and insula- 
tor states and should be described by itinerant-electron 
models rather than by localised-electron ones. Since the 
metal-insulator transition in the model under considera- 
tion turns out to be strongly first-order type, the gap is 
not small. The value of the gap depends strongly on the 
inter-site Coulomb repulsion parameter V, as is demon- 
strated in Fig. 4. The insulating properties cannot be 
described correctly neglecting the charge fluctuations on 
the rung (too small values of V or using single-site DMFT 
approach instead of the cluster one). 

To conclude, we present for the first time results on 
the cluster DMFT calculations for the insulating phase 
of NaV2 05 above charge-ordering transition using the 
extended Hubbard model with the first-principle tight- 
binding parameters. The non-local charge fluctuations 
and inter-site Coulomb interaction is of crucial impor- 
tance for the formation of insulating state. 
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